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Abstract. The method developed by Vinti for predicting the motion of an 
artificial satellite is used to obtain the first order partial derivatives in the 
conditional equations. These in turn are used in orbit improvement pro- 
grams to determine the mean elements of Vinti's theory, corresponding 
to a separable Hamiltonian. The results have been programmed and 
tested on the IBM 7094 Mod. II. 

INTRODUCTION 

The kinetic equations of motion (Vinti 1959) constitute a reduction of 
the problem of artificial satellite motion to quadratures. Solution for the 
reference orbit consists of inverting these equations to find the coordi- 
nates p and 7 } as functions of time. The <x. and (i = 1,2,3) are the 
Jacobi constants. If the initial conditions are known, one can readily 
evaluate the a.. Then if one can evaluate the integrals, one can evaluate the 
/3. . Evaluating the integrals and consequently performing the required inver- 
sion depend on factoring the quartics F(p) and G (rj) . These are given by 

F(p) = c 2 a 2 + (p 2 + c 2 ) (-a 2 + 2 pp + 2 a, p 2 ) 

- (- 2 ap (p - pj) (p 2 - p) (p 2 + A p + B) 

0 ( 77 ) r. - a 2 + (1 - ij 2 ) (a 2 + 2 a, c 2 rj 2 ) 

- - 2 a, c 2 (t 7 2 _ rj 2 ) (tj 2 - rj 2 ) 

Comparison of the two expressions for F(p) leads to four algebraic 
equations for the quantities Pj , p 2 , A , and B in terms of a j, a 2 , a 3 » ^ » 
and c 2 . Here we use the constants of the motions f\ , > P$ • an( ^ 
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obtainable from initial conditions. One can then find r/ 0 and rj 2 explicitly, 
and p , p , A , and B, by successive approximations through order J 2 , 
where J 2 is the coefficient of the second zonal harmonic of the planet's 
gravitational potential. 


There is another set of elements introduced by Izsak (Izsak I960) 
which result in the immediate factoring of the quartics F (p) and G(tj) • 
These are /3j , ft 2 , /? 3 , and 

a = \ (Pi +/° 2 ) * 

( P 2 - P.) 

(p 2 + p.) ’ 

I = sin -1 v Q 

It is possible to express a x , a 2 , a 3 , n 0 , v 2 , p t , p 2 , A, and B explicitly 
in terms of a, e , and I. To obtain a , e, and I, we first approximate 
them by a Q , e 0 , and i Q , compute the orbit by the above method, compare 
results with observations over an interval of time, find residuals, and 
then do iterated differential corrections (Bonavito 1962, 1964). 
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Since the integrands contain double valued functions, Vinti has intro- 
duced a set of uniformizing variables E, v, and t p, to simplify them. 

These are defined as variables analogous to the eccentric anomaly, true 
anomaly, and argument of latitude respectively. 

Let L Q denote any observed quantity characteristic of the method of 
tracking of the satellite. For example, this can represent a direction 
cosine for the minitrack radio system, a range rate for doppler tracking, 
or a right ascension or declination for optical tracking. Denoting a, e, 
v 0 , /3j , JS 2 , and fi i by q. (i = 1 , 2, . . . 6), we can then write the conditional 
equations as 



Here L c represents the computed value of the observed quantity. The 
first derivatives are given by the equation 

dL dL 3x 3L By dL Bz 

C cm r J m cm 

Bn dx da. Bv 3a. d z 3q. 

With use of the appropriate transformations, the expressions obtained 
for L , BL /3x , dL /By , and 3L / <)z , will depend only on the local 

coordinates. One can then determine them from the values of the inertial 
coordinates x , y, z, computed by the Vinti orbit generator for specific 
observation times. 
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By means of the equations 

P = &( 1 - e cos E), 

V - V 0 s i n 

the inertial coordinates can be written 

x j/[a 2 (1 - e cos E ) 2 + c 2 ] (1 - 17 2 Sin 2 < p) cos <£, 
y - y/(a 2 (1 - e cos E ) 2 ■<- c 2 ] (1 - p 2 sin 2 \p) sin <p, 

and 

z = a(l - ecosE) 17 sin <//. 

One can derive the formulas for 3x m /dq. , 3y m /3qj, and dz m /Bq. in 
terms of p, r/ , and <p by direct differentiation of x, y , z , with respect 
to a, e, i / 0 , fi j , , and /' 3 . Since the coordinates x , y , and z are 

functions of the uniformizing variables, it is necessary to obtain the 
expressions for dE/dq. , and di/i/'dq. (Bonavito 1964). 

This can be performed in the following manner: 

Differentiating the expressions on page 199 of reference 5 (Vinti 1961), 
we list 

1 . _2 and — — 
dt »i 

aE o 

from the Kepler equation. 

dq. 
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3 # by using v = M + v 0 , E - + E fl , and the anomaly connec- 

<)q 


4. 


tions. 

V 0 ,)M i 


, and 


BE, 


Bq, Bq. Bq. 


Bv. 


5. — l by using v = M + v„ + v. , E = M f + E 0 + Ej , and the anomaly 

3q. 


connections. 


Bi/ij c3M 2 3E 2 9E B(M s + E 0 + E x + E 2 ) 

6. - — , - — , and ^ — , which gives =r — = 

3q. ' Bq.’ Bq.’ g Bq. 
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Bv, 


7. Then — — - with use of v - + v 0 + v j + v 2 , E = M s + E Q + 

dq, 


Ej + E 2 , and the anomaly connections, which gives 
Bv + v o + v i + v 2 ^ 


Bq. 


Bq ; 


B'/', 


8. Then — — yields 

Bq ; Bq. 


_ a V, + >/'o + V'- 1 ! + V 2 > 


Bq x 


9. ,^X , where x is an angle that equals p whenever is a multiple 

Bq; 


of — , and which also satisfies the equation 


tan x = (1 - V) 1 2 tan ^ 


10* Finally 


Bq ; 
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Sinn' the quantities JE/<)q Oi/>/<)q. ( and *>./>/ >)q. are obtained only from 
the theoretical calculation, one can then immediately determine the first 

derivatives for use with any method of tracking whatsoever after deciding j 

on the form of L and its partial derivative with respect to the local 
coordinates (Bonavito 1964). Corrections to the orbital elements Aq. 
are computed by an iterated least square procedure after a sufficient 
number of observations are obtained, a minimum of six conditional equa- 
tions being necessary to determine the problem. Classically speaking, 
the initial coordinates of the system point are shifted in p, q phase 
space such that the continuous evolution of a canonical transformation 
with this new value of the Hamiltonian produces a motion of a mechanical 
system consistent with observation, in this case the satellite orbit. 

» 

In the following, cr = i/a(D 0 -D c ) 2 /'(N 1 - N 2 ) where D Q and D c are the j 

observed and computed comparison parameters (i.e.: direction cosines, 

right ascension, declination, etc.), Nj is the number of equations of con- 
dition (20), and N 2 is the number of variables (6). Figure 1 is a plot over 
an allowed ten iterations of the standard deviation of fit a, times 10 3 for 
direction cosine observation data from the Relay II Satellite (#64031). 

The criterion for accepting observations is AL (AM) < 3 X 10' 3 . The 
twenty equations of condition (without weighting factors) used in a least 
squares routine to find the six Izsak elements, extend over an arc of 


three hours. Figure 2 is a similar plot for Smithsonian Astrophysical 
Observatory optical data from the ANNA Satellite. Twenty equations of 
condition are used for an arc of seventy five hours following injection, 
and the criterion for accepting observations is AL (AM) < 20 X 10‘ 3 
milliradians. Although the number of iterations was extended to 
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ten, the program converged after the third iteration for the Relay 
II data, and after the second iteration for the data from Satellite 
ANNA. The value of the eccentricity for each of the above satellites 
is considerably different. Thus, for the initial epoch, they were 
taken to be 0.23597617 for Relay II, and 0.00671710 for ANNA. 

Figure 3 is a plot for the first sever, equations of condition of Satellite 
ANNA covering a period of the first forty four hours following injection. 
While a convergence criterion of 0.2 x 10~ 3 for the standard deviation of 
fit was met using the ANNA data and twenty equations of condition, the 
theory was able to meet a value of 0-04 x 10~ for the case of seven 
equations of condition. Also, it has been shown that the Vinti orbit 
generating program computes positions and velocities very rapidly and 
with great accuracy (Bonavito 1962). For seventy equations of condition 
of the Relay II Satellite extending over a five hour period, the program 
converged on the third iteration to a standard deviation of fit criterion of 
2.7 x 10 within thirty seconds. All ol the above investigations were 
conducted on (he IBM 7094 electronic digital computer. These 1 * results 
are based only on the Vinti reference: potential, without the effects of the 
third harmonic or the residual fourth harmonic. Inclusion of these other 
effects, which could be accomplished by use of the perturbation theory 
of Vinti (1963) would be expected to furnish computed orbits which woulei 
agree with observation over a longer interval of time. 
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